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Local-density-approximation calculations are used to show that the metal-metal 
distance along the chains controls the charge-density wave (CDW) in halogen-bridged 
transition-metal linear-chain (MX) compounds. The strength of the CDW can be 
understood in terms of a two band Su-Schrieffer-Heeger model if a hard-core ion-ion 
repulsion potential is also added. We predict a second-order phase transition from an 
insulating to a semi-metallic ground state and explain trends in dimerization, bond- 
length ratios, band gaps, and Raman breathing modes in terms of the metal-metal 
distance. 
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During the last few years halogen-bridged transition-metal linear-chain (MX) compounds 
have been extensively studied 0,0,0,^J^J^,0,|],^|,|1O|,|TT],|T2|,|T3[. One reason for the increased 
interest in these materials is that their various competing ground-states such as charge- 
density-wave and spin- density- wave are sensitive to tuning by chemical substitutions, pres- 
sure, or doping P|,|2]J^,^,p||5|,[7||8|,|9|, pi]|jn| . In particular, progress has been achieved recently 
and independently by groups in Japan anc ^ m Los Alamos |9|, ^0| , p]| in tuning the 



charge density wave (CDW) strength over a wide range by using different counter-ions and 
equatorial ligands. 

In this Letter, we use local-density-approximation (LDA) full-potential linear muffm-tin- 
orbital (FPLMTO) calculations to study the tuning of the CDW of the MX systems as a 
function of the metal-metal (M-M) distance along the chain. By comparing experimental 
results for MX (X= Br, CI) systems with varying M-M distances to our calculated results 
for Pt2Xe(NH3)4 (X= Br, CI) under uniaxial stress, we show that the key parameter is the 
M-M distance. Thus, the effect of the different ligands on the properties of these systems 
is indirect; they tune the M-M distance which is really responsible for the change in the 
electronic properties. Moreover, we also demonstrate that an extended Su-Schrieffer-Heeger 
(SSH) model with a universal set of parameters fits the entire class of Pt charge-density-wave 
MX materials. Thus, we have now provided a crucial theoretical underpinning that explains 
the CDW MX systems. 

Our principal results for the uniaxial stress are: (i) The dimerization, the condensation 
energy, the Raman breathing mode, and the band gap are gradually reduced, and vanish 
when the Pt-Pt distance is reduced by nearly 7% from its ground-state value, (ii) Under 
the same conditions, we predict that the PtBr MX system exhibits a second-order phase 
transition from an insulator to a semi-metal, (iii) The two-band SSH model |TJ[ describes 
the tuning of the CDW strength only if we also introduce a repulsive potential between the 
metal and halogen ions. 

The calculations presented in this paper are based on the LDA with an all-electron 
scalar-relativistic FPLMTO basis-set as described in our previous work |y|. This method 



was successfully used to show that the dimerization and the insulating ground state of 



Pt 2 Br 6 (NH 3 )4 are correctly predicted provided that we include the full ligand structure [|T3 
Under uniaxial stress, we rigidly move the ligands with the metal atoms. 

A single MX chain (M = Ni, Pd, or Pt; and X = CI, Br, or I) consists of alternating MX 
units, arranged such that two neighboring X atoms move closer to one of the M atoms [|T],3 . 
This dimerization is schematically represented by ...M 3 ~ 5 ...XM 3+<5 X...M 3 ~ <5 ...XM 3+,5 X... Due 
to the strong electronegativity of the X atoms, the dimerization causes an alternation of the 
valence character of the metal atoms between 3 — 5 and 3 + 5. The strength of the CDW 
produced is then directly proportional to the dimerization. 

Metal-insulator phase transition. Figure [I] shows the LDA total energy of the PtBr system 
for Pt-Pt distances less than the experimental lattice parameter (5.55 A) as a function of 
the dimerization £. Under uniaxial stress, the double-well structure in the energy collapses 
to a single well, and the condensation energy (the difference between the energy at zero 
dimerization and the minimum energy) and the dimerization decreases as the Pt-Pt distance 
is reduced. This agrees with earlier work by Sakai et al. on Pt(en)Cl(C104)2 under pressure 
||. At a Pt-Pt distance of 5.18 A (which is equivalent to a uniaxial stress of 89 kbars based 
on the first derivative of the total energy versus cell volume) the condensation energy and the 
dimerization vanish simultaneously, and the system exhibits a second-order phase transition 
from an insulating to a semi-metallic ground-state. This is confirmed by the calculated 
band gap which decreases gradually and closes when the Pt-Pt distance is near 5.18 A (see 
the right side of Fig. @). The possibility of this transition in [Pt(NH 3 )4XPt(NH 3 )4X] 4+ 
chains was first suggested by Whangbo and Foshee using an extended Huckel method [[|. 
It is also consistent with pressure measurements on the PdBr system, where 140 kbar of 
hydrostatic pressure caused a nine-order-of- magnitude increase in electrical conductivity || , 
due to an increased carrier concentration (because of the reduced gap). However, beyond 
140 kbar of pressure a leveling off of the conductivity versus pressure and an increase of 
the activation energy was observed upon application of pressure up to 350 kbar ||. Similar 
conclusions were reached by other researchers [HHH for Pd(chxn)Br(Br)2, Pt(en)Cl(C104)2, 



and Pt(en)Br(C104)2 systems respectively. 

Following the suggestion of Interrante et al. || that a compression of the metal-ligand 
bond length may increase the band-gap, we have performed self-consistent calculations of the 
PtCl band structure in which the interchain distance is decreased by 5% without modification 
of the M-M distance of 5.45 A along the chain. When the ligand-metal distance is held 
constant the band gap is unaffected, whereas a 5% decrease in the ligand-metal bond-length 
increases the energy band gap by 0.25 eV. These results coupled with the decrease in the 
band gap when the metal-metal distance along the chain decreases suggest an intricate 
interplay between the ligand-metal distance and the conduction process along the chain, 
and may serve as a qualitative explanation for the increase of the band-gap beyond some 
critical pressure. Thus, real pressure is a more complicated situation than the uniaxial 
pressure generated by chemical subsitutions, where the ligand-metal distance is constant, 
the CDW state is controlled by the M-M distance, and a metallic state should occur as 
shown by our calculation. 

Softening of the Raman mode. The use of a harmonic oscillator model in the effec- 
tive potential given by the total energy curve allows us to estimate the optical (Raman) 
breathing-mode phonon frequency as a function of the Pt-Pt distance. This breathing mode 
is an oscillation of the two Br chain atoms around a Pt atom at rest. Since the curvature in 
the total energy curve (at its minimum) decreases under uniaxial pressure, this frequency is 
expected to decrease. 

Figure displays our calculated Raman phonon frequencies and band g function 
of the Pt-Pt distance. We observe an almost linear softening of the Raman mode as the 
M-M distance is shortened. In the vicinity of the phase-transition, where the curvature 
of the total energy becomes flat (cf. Fig. [1]) the phonon mode softens considerably. The 
phonon frequencies for Pt-Pt distances below 5.37 A are not displayed, because the use of 
a simple harmonic oscillator to calculate the phonons is dubious since the lattice distortion 
is comparable to the zero-point-motion, and the phonon frequency is comparable to the 
condensation energy. In this case, thermal effects and zero-point motion, which may cause 
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a breakdown in the Born-Oppenheimer approximation, must be taken into account. 

Figure |3] gives a comparison between the calculated short Pt 3+<5 -X and long Pt 3_<5 -X 
bond lengths [X = Br (open circles), CI (filled circles)] versus the Pt-Pt lattice parameter 
with the different experimental results obtained for different MBr/Cl compounds (M = Pt, 
Pd, and Ni). The different ligands and counter- ions produce a range of M-M distances for 
the different compounds. For example, experimentally, substitution of the CIO4 counter-ion 
by Br decreases the metal-to-metal distance along the chain because the hydrogen bond 
between NH2 of ethylenediamine (en) or (lR,2R)-cyclohexanediamine (chxn) ligands and 
Br counter-ion is strengthened J7| . The substitution of (chxn) ligand by (en) ligand reduces 
the metal-metal distance because (chxn) has a non-planar structure whereas (en) is planar 
i- 

The experimental results of all these different materials track the calculations under uni- 
axial pressure. This strongly suggests that the M-M distance is the key parameter controlling 
the CDW strength. 

In addition, the short M 3+<5 -X bond distance remains essentially constant throughout 
the series of compounds presented, revealing the contraction of the long M 3_<5 -X bond as 
responsible for the large change in the dimerization, therefore reducing the strength of the 
CDW. Our calculations describe the experimental situation rather well, even though the 
experimental results are mostly for different MX chain systems with different ligands. But 
this surprising agreement can be understood qualitatively as follows. In our previous work, 
we have shown that the dimerization in PtBr system is driven by an electron-phonon coupling 



confined to a single chain |T3[ . This strong electron-phonon coupling reduces the short M-X 
distance to the sum of the ionic radii of the metal and the bridging halogen ions [[/J. Thus, 
the effect of uniaxial stress is to decrease the long M-X bond length which is much softer 
than the shorter one. To reduce the short bond length, one has to act against the strongly 
anharmonic potentials of the M 3+<s and X ions Jl6|] . 

We have analyzed quantitatively the effects of uniaxial stress on the strength of the CDW 



in terms of a two-band single-chain SSH model |12 



#o = £{Ho + oA,)(cU +lff + ^-c.) + [(-l)'eo - A (A, + A^)]^ J 

+ E 4?A? 1+i) + E «^(A a + A a+1 )^ , (1) 

1,3=0,1 i,J=0,2 

where cj CT (q j(J ) creates (annihilates) an electron at site x\ with spin a, and where M and 
X occupy even and odd sites, respectively. The electron-phonon parameters a, foi = Pm, 
and foi+i = fix couple the lattice (through A ; = — x{) to the electronic structure. The 
energies ±eo and the to are the M and X energy levels and the hopping integral at zero 
dimerization, respectively. The spring constants Kyjx between M and X and Kyj M between 
M and M model represent elastic energies that model the rest of the electronic structure 
that is not included explicitly in the one-dimensional electronic structure |l9fl . 

The bands in the vicinity of the Fermi energy deviate slightly from the nearest-neighbor 
tight-binding bands of the SSH model, Eq. (1). Nonetheless, this model, with only one set 
of parameters, fits nicely the relevant bands for any Pt-Pt distance. 

Anharmonic potential. One failure of the SSH model is its prediction of an almost con- 
stant ratio of the dimerization to the M-M distance (about 3.6% for PtBr) for all M-M 
distances, in disagreement with experiment and LDA (which predicts a vanishing dimer- 
ization for a sufficiently short M-M distance). We have extended the model by adding a 
hard-core ion-ion repulsive potential between the metal and halogen atoms along the chain. 
For simplicity, we have used a next-neighbor repulsive potential. This changes the model 
Hamiltonian to 
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(itfe + a, - R c y 



where A and q are constants representing the strength and hardness of the potential, Rj$x 
is the M-X distance at the experimental equilibrium ground-state at zero dimerization, and 
R c is a hard-core distance between M and X [p| . 

The Hamiltonian of Eq. (fj) fits the LDA total energy of PtBr under uniaxial stress rather 
well (see Fig. [TJ, dashed curves). The parameters used to describe the LDA total energies are 
|r?| given in Table |I[ In Fig. ^| we have also depicted the short M 3+<5 -X and the long M 3_5 -X 



bond lengths predicted by the model given by Eq. (2) as a function of the Pt-Pt distance. 
Notice that the model is in good agreement both with the LDA and the experimental results. 
So, the constant M 3+,5 -X short bond is explained as a result of the competition between a 
strong electron-phonon coupling acting along the chain and a hard-core repulsive potential 
between the metal and halogen ions. 

To further illustrate the model given by Eq. (2), we have calculated the short and long 
bond for PtCl. The LDA is used to get the dimerization for the ground state (Pt-Pt distance 
of 5.27 A ||20||)) then the model is fitted to the LDA results with the requirement that 
the short bond-length is equal to the LDA prediction. This procedure provided the model 
parameters shown in Table |[ Figure |] shows that our prediction for PtCl compares favorably 
with the available experimental results. The model also predicts a metallization for PtCl at 
a Pt-Pt distance of 4.78 A. Notice that the Ni-X sytems (numbers 6 and 13 in the figure) 
don't dimerize; in fact, they are magnetic, which opens a gap that prevents them from being 
metallic. 

In conclusion, we have found that LDA correctly predicts the variation in CDW strength 
of PtBr and PtCl under uniaxial stress, and that the key parameter is the M-M distance. 
We have shown that the dimerization, the band gap, and the breathing phonon mode all 
decrease linearly with the M-M distance in agreement with experiments for related systems 
PJ5|fl;HflR TOPT| ■ At a Pt-Pt distance close to the sum of the ionic radii a second-order 
phase transition to a metallic state appears. 

Finally, we have shown that the SSH model required the addition of a hard-core repulsive 
potential between the metal and the halogen ions to reproduce the correct tuning of the 
CDW strength in the MX systems. This shows that the constant short M 3+<s -X distance is 
an equilibrium distance resulting from a competing strong electron-phonon coupling and a 
hard-core repulsive potential between the metal and the halogen ions. 

We thank R. H. Mckenzie for stimulating discussion. Partial support was provided by 
the Department of Energy (DOE) - Basic Energy Sciences, Division of Materials Sciences. 
Supercomputer time was provided by the Ohio Supercomputer and by the DOE. 
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FIGURES 

FIG. 1. LDA total energy of PtBr as a function of dimerization ratio £ (the ratio of half the 
difference between the two M-X distances to the M-M distance in %) for different Pt-Pt distances. 
The Pt-Pt distance are, from the upper curve to the lower one, 5.18 A, 5.29 A, 5.37 A, 5.45 A, and 
5.55 A. The filled circles are the calculated LDA total energies and the full curves are obtained 
using a polynomial least-squares fit of order four. The dashed curves represent the total energies as 
obtained from the Hamiltonian given by Eq. (2) (see text). An energy offset of 0.05 eV instead of 
the calculated energy offsets has been added between the different curves for clarity of presentation 



[15]. Under uniaxial stress the size of the dimerization and the dimerization energy are reduced 
and, at Pt-Pt distance of 5.18 A the dimerization vanishes and a second-order phase transition to 
a metallic state appears. 

FIG. 2. Calculated (open circles) and measured (l[ (filled circle) Raman breathing mode, and 
calculated band gap energy (filled triangles) as a function of the Pt-Pt distance (the experimental 
band gap under pressure decreases [jll| but the absolute magnitude is not yet known). The linear 
softening of the breathing phonon mode is in good agreement with the experimetal results of Ref. 
[§.11] and the results (long dash) of the model given by Eq. (2). 



FIG. 3. LDA calculated short M 3+,5 -X and long M 3 ~ 5 -X bond lengths for Pt 2 Br 6 (NH 3 )4 
(open circles) and Pt 2 Clg(NH3) 4 (filled circles) as a function of M-M distance. The filled tri- 
angle symbols represent experimental data for various MX systems |JT|: l-Pt 2 Br 6 (NH 3 ) 4 ; 
2-[Pt 2 (en) 4 Br 2 ] (C10 4 ) 4 ; 3-[Pd 2 (en) 4 Br 2 ] (C10 4 ) 4 ; 4-[Pt 2 (chxn) 4 Br 2 ] (Br) 4 ; 5-[Pd 2 (chxn) 4 Br 2 ] 
(Br) 4 ; 6-[Ni 2 (chxn) 4 Br 2 ] (Br) 4 ; 7-[Pt 2 (chxn) 4 Cl 2 ] (C10 4 ) 4 ; 8-[Pt 2 Cl 2 (NH 3 ) 8 ] (HS0 4 ) 4 ; 
9-[Pt 2 (en) 4 Cl 2 ] (C10 4 ) 4 ; 10-[Pt 2 (tn) 4 Cl 2 ] (BF 4 ) 4 ; ll-[Pt 2 (en) 4 Cl 2 ] (Cu0 4 ) 4 ; 12-[Pt 2 (chxn) 4 Cl 2 ] 
(Cl) 4 ; 13-[Ni 2 (chxn) 4 Cl 2 ] (Cl) 4 . The solid and dashed curves are guides to the eye. The dot-dashed 
curve is the results predicted by the model of Eq. 2. Notice that the calculated and measured short 
M-X bond lengths are roughly independent of the M-M distance and that the agreement between 
the calculations and the experimental results is good. 
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TABLES 

TABLE I. Parameters of the model given by Eq. 2 obtained from the fit to the LDA results. 
The units of 2eo and to is eV, that of a is eV/A, and that of is eV/A^ 2 . The j3 parameters 
are found to be zero. 
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